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Abstract. The symmetry energy is an important quantity in the equation of 
state of isospin asymmetric nuclear matter. This currently unknown quantity 
is key to understanding the structure of systems as diverse as the neutron- 
rich nuclei and neutron stars. At TAMU, we have carried out studies, aimed 
at understanding the symmetry energy, in a variety of reactions such as, the 
multifragmentation of 40 Ar, 40 Ca + 58 Fe, 58 Ni and 58 Ni, 58 Fe + 58 Ni, 58 Fe reac- 
tions at 25 - 53 AMeV, and deep-inelastic reactions of 86 Kr + 124 ' 112 Sn, 64 > 58 Ni 
(25 AMeV), 64 Ni + 64 > 58 Ni, 112 > 124 Sn, 232 Th, 208 Pb (25 AMeV) and 136 Xe 
+ 64 ' 58 Ni, 112 < 124 Sn, 232 Th, 197 Au (20 AMeV). Here we present an overview 
of some of the results obtained from these studies. The results are analyzed 
within the framework of statistical and dynamical models, and have important 
implications for future experiments using beams of neutron-rich nuclei. 
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1. Introduction 

The equation of state of isospin asymmetric nuclear matter is a fundamental quan- 
tity that determines the properties of systems as small and light as an atomic 
nucleus, and as large and heavy as a neutron star. The key unknown in the EOS 
of asymmetric nuclear matter is the symmetry energy. Recently the possibility of 
extracting information on the symmetry energy and the isospin (neutron-to-proton 
ratio) of the fragments in a multifragmentation reaction has gained tremendeous 
importance 12]- Such information is of importance for understanding key prob- 
lems in astrophysics [23, and various aspects of nuclear physics such as the structure 
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of exotic nuclei (the binding energy and rms radii) [|IJ |3] and the dynamics of heavy 
ion collisions 

Traditionally, the symmetry energy of nuclei has been extracted by fitting the 
binding energy of the ground state with various versions of the liquid drop mass 
formula. The properties of nuclear matter are then determined by theoretically 
extrapolating the nuclear models designed to study the structure of real nuclei. 
However, real nuclei are cold, nearly symmetric (N « Z) and found at equilibrium 
density. It is not known how the symmetry energy behaves at temperatures, isospin 
(neutron-to-proton ratio) and densities away from the normal nuclear matter. The- 
oretical many-body calculations [ 7 and those from the empirical liquid drop mass 
formula predict symmetry energy near normal nuclear density (« 0.17 fm -3 ) 
and temperature (T w MeV), to be around 28 - 32 MeV. 

In a multifragmentation reaction, an excited nucleus expands to a sub-nuclear 
density and disintegrates into various light and heavy fragments. The fragments are 
highly excited and neutron-rich ; their yields depend on the available free energy, 
which in turn depends on the strength of the symmetry energy and the extent to 
which the fragments expand. By studying the isotopic yield distribution of these 
fragments, one can extract important information about the symmetry energy and 
the properties of the fragments at densities, excitation energies and isospin away 
from those of ground state nuclei. Here we present some of the results obtained from 
various measurements carried out at the Cyclotron Institute of Texas A&M Uni- 
versity (TAMU). We present these results in the framework of both, the statistical 
and the dynamical multifragmentation models. 

2. Isoscaling and Symmetry Energy 

In a multifragmentation reaction, the ratio of isotope yields in two different re- 
actions, 1 and 2, R2i(N,Z) = Y 2 (N, Z)/Yi(N, Z), has been shown to obey an 
exponential dependence on the neutron number (TV) and the proton number (Z) of 
the isotopes, an observation known as isoscaling ^3 The dependence is 
characterized by a simple relation, 



where, Y 2 and Y\ are the yields from the neutron-rich and neutron-deficient 
systems, respectively. C is an overall normalization factor, and a and f3 are the 
parameters characterizing the isoscaling behavior. 

Theoretically, isoscaling has been predicted by both, statistical [E3] as well as 
dynamical [Ej multifragmentation models. In these models, the difference in the 
chemical potential of systems with different N/Z is directly related to the scaling 
parameter a. The scaling parameter a is proportional to the symmetry energy 
through the relation, 



R 2 i(N, Z) = Y 2 {N, Z)/Yi(N, Z) = Cexp{aN + (3Z) 



(1) 




(2) 
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Fig. 1. Isotopic yield ratios for 40 Ar, 40 Ca + 58 Fe, 58 Ni reactions at 25 - 53 
MeV/nucleon, and deep-inelastic reactions of 64 Ni + 64 - 58 Ni, 112 < 124 Sn, 232 Th, 208 Pb 
(25 MeV/nucleon) 



In the above equation, Z\ , Ai and Z 2 , A 2 are the charge and the mass numbers 
of the fragmenting systems, T is the temperature of the system and C sym , the 
symmetry energy. The parameter a, has been shown to be independent of the 
complex nature of the secondary de-excitation of the primary fragments, and is 
thus a robust observable for studying the symmetry energy [110]. While it is well 
established that many versions of statistical models show very little or no difference 
between the a values for the primary and the secondary fragments, the same may 
not be true for the a values obtained from dynamical models. The origin of this 
discrepancy between the two approaches is currently being debated and not fully 
understood. It will be shown from the present study that the difference between the 
primary and the secondary en's in statistical model is very small. We will presume 
the difference to be insignificant in the dynamical model framework. 

Fig. 1 shows the experimentally determined isotopic yield ratio as a function 
of neutron number N, for some of the reactions studied using beams from the 
K500 Cyclotron at TAMU. The figure on the left shows the ratios for the ( 40 Ar + 
58 Fe)/( 40 Ca + 58 Ni) and ( 40 Ar + 58 Ni)/( 40 Ca + 58 Ni) pairs of reactions. The one 
on the right is for the deep-inelastic reactions of 64 Ni + Ni, Sn, Th, Pb at 20 AMeV. 
One observes that the ratios for various elements in a given reaction pair lie along 
a straight line in the logarithmic plot and align with the neighboring elements in 
accordance with the relation given in equation 1. This feature is observed for all 
the beam energies and the pairs of reactions studied. The alignment of the data 
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points varies with beam energies as well as the pairs of reaction. The ratio for each 
elements (Z) were simultaneously fitted using an exponential relation (shown by the 
solid lines) to obtain the slope parameter a. In the following sections, we use these 
experimentally determined a's to study the symmetry energy using the statistical 
and dynamical model interpretation of the multifragmentation reaction. 

3. Symmetry Energy from a Statistical Model Approach 

The Statistical Multifragmentation Model (SMM) [E3 d is the most widely used 
model for describing multfragmentation reactions. It is based on the assumption 
of statistical equilibrium at a low density freeze-out stage. All breakup channels 
composed of nucleons and excited fragments are taken into account and considered 
as partitions. During each partition the conservation of mass, charge, energy and 
angular momentum is taken into account, and the partitions are sampled uniformly 
in the phase space according to their statistical weights using Monte Carlo sampling. 
The Coulomb interaction between the fragments is treated in the Wigner-Seitz ap- 
proximation. Light fragments with mass number A < 4 are considered as elementary 
particles with only translational degrees of freedom ( "nuclear gas" ) . Fragments with 
A > 4 are treated as heated nuclear liquid drops, and their individual free energies 
Fa,z are parametrized as a sum of the volume, surface, Coulomb and symmetry 
energy, 

Fa,z = F\ z + F% z + E AZ + E S A V ™ (3) 

where F AZ = (—W — T 2 /e )A, with parameter e related to the level density 
and W = 16 MeV being the binding energy of infinite nuclear matter. F A z = 
B A 2 ^ 3 [(T 2 - T 2 )/{T 2 + T 2 )] 5 / 4 , with B = 18 MeV being the surface co-efficient 
and T c = 18 MeV being the critical temperature of infinite nuclear matter. E AZ = 

cZ 2 /A 1/3 , where c = (3/5)(e 2 /r c )[l - (p/p ) 1/3 } is the Coulomb parameter ob- 
tained in the Wigner-Seitz approximation with charge unit e, and r Q — 1.17 fm. 
e az = C sym (A - 2Z) 2 /A, where C s ym — 25 MeV is the symmetry energy co- 
efficient. These parameters are adopted from the Bethe-Weizsacker mass formula 
and correspond to the assumption of isolated fragments with normal density in the 
freeze-out configuration. The value of the symmetry energy co-efficient C sym is 
taken from the fit to the binding energies of isolated cold nuclei in their ground 
states. In a multifragmentation process the primary fragments are not only excited 
but also expanded. 

Figure 2 shows a comparison between the SMM calculated and the experimen- 
tally observed values of a. The left side of the figure corresponds to the ( 40 Ar + 
58 Ni)/( 40 Ca + 58 Ni) and the ( 40 Ar + 58 Fe)/( 40 Ca + 58 Ni) pairs of reactions. The 
one on the right corresponds to ( 58 Fe + 58 Fe)/( 58 Ni + 58 Ni) and ( 58 Fe + 58 Ni)/( 58 Ni 
+ 58 Ni) pairs of reactions. The dotted lines in the a versus excitation energy plot 
corresponds to a calculated from the primary fragment distribution and the solid 
lines to those calculated from the secondary fragment distribution. The symbols 
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Fig. 2. Statistical model comparison of the isoscaling parameter with the experi- 
mentally determined values as a function of the excitation energy. 



correspond to the experimentally determined a's. The figure on the left clearly 
shows that the experimentally determined a's are significantly lower than the cal- 
culated values using the standard value of the symmetry energy, C sym = 25 MeV. 
To explain the observed dependence of the isoscaling parameter a on excitation 
energy, the C sym of the hot primary fragment in the SMM calculation was varied 
in the range 25 - 15 MeV. As shown in the center and the bottom panel of the left 
figure, the isoscaling parameter decreases slowly with decreasing symmetry energy. 
The experimentally determined a can be reproduced for both pairs of systems at all 
excitation energies using a symmetry energy value of C sym — 15 MeV. This value of 
the symmetry energy is significantly lower than the value of C sy m — 25 MeV often 
used for the isolated cold nuclei in their ground states. On the right side of the 
figure, we show the comparisons for the ( 58 Fe + 58 Fe)/( 58 Ni + 58 Ni) and ( 58 Fe + 
58 Ni)/( 58 Ni + 58 Ni) pairs of reactions. Once again a lower value of symmetry energy 
C syrn is required to explain the experimental data. Furthermore, one also observes 
a small dependence of the symmetry energy with increasing excitation energy. 

A similar behavior of the C sym is also observed from the deep-inelastic reactions 
studies of 86 Kr + 124 < 112 Sn, 64 > 58 Ni (25 AMeV), 64 Ni + Ni, Sn, Th, Pb (25 AMeV) 
and 136 Xe + Ni, Sn, Th, Au (20 AMeV) [|T3j. Fig. 3 shows the values of the C sym 
obtained from these reactions. The values were obtained using equation 2 with two 
different assumptions for temperature determination ; the Fermi gas temperature 
(closed symbols) and the expanding mononucleus temperatures (open symbols). The 
symmetry values for both set of temperatures show decreasing trend with increasing 
excitation energy. 

The above comparison of the experimentally observed isoscaling parameter with 
the statistical multifragmentation model therefore shows that a significantly lower 
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Fig. 3. Symmetry energy as a function of the excitation energy for 86 Kr + 124,112 Sn 
(circles), 64 > 58 Ni (squares) (25 AMeV), 64 Ni + Ni, Sn, Th, Pb (triangles) (25 AMeV) 
and 136 Xe + Ni, Sn, Th, Au (inverted triangles) (20 AMeV) reactions. 

value of the symmetry energy is required to explain the isotopic composition of the 
fragments produced in a fragmentation reaction. This indicates that the properties 
of nuclei at high excitation energy, isospin and reduced density are very sensitive to 
the symmetry energy. Similar hot and neutron-rich nuclei are routinely produced in 
the interior of a collapsing star and subsequent supernova explosion, where a slight 
decrease in the symmetry energy can significantly alter the elemental abundance 
and the synthesis of heavy elements [ llri| . The present observations can provide 
important inputs for the understanding of the nuclear composition of supernova 
matter. 

4. Symmetry Energy from a Dynamical Model Approach 

In the following, we analyze the above results in the dynamical model framework 
using the Anti-symmetrized Molecular Dynamic (AMD) calculation AMD 
is a microscopic model that simulates the time evolution of a nuclear collision, 
where the colliding system is represented in terms of a fully antisymmetrized prod- 
uct of Gaussian wave packets. During the evolution, the wave packet centroids 
move according to the deterministic equation of motion. The followed state of the 
simulation branches stochastically and successively into a huge number of reaction 
channels. The interactions are parameterized in terms of an effective force acting 
between nucleons and the nuclcon-nuclcon collision cross-sections. The beauty of 
the dynamical models is that it allows one to understand the functional form of the 
density dependence of the symmetry energy at a very fundamental level i.e., from 
the basic nucleon-nucleon interactions. Theoretical studies [[7] based on microscopic 
many-body calculations and phenomenological approaches predict various forms of 
the density dependence of the symmetry energy. In general, two different forms 
have been identified. One, where the symmetry energy increases monotonically 
with increasing density (" stiff " dependence) and the other, where the symmetry 
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energy increases initially up to normal nuclear density and then decreases at higher 
densities (" soft " dependence). 

Determining the exact form of the density dependence of the symmetry energy 
is important for studying the structure of neutron-rich nuclei [EI 111 j an( i studies 
relevant to astrophysical problems, such as the structure of neutron stars and the 
dynamics of supernova collapse [Ej- For example, a " stiff " density dependence of 
the symmetry energy is predicted to lead to a large neutron skin thickness compared 
to a " soft " dependence [EI- Similarly, a " stiff " dependence of the symmetry 
energy can result in rapid cooling of a neutron star, and a larger neutron star radius, 
compared to a soft density dependence [H5|- 

Recently, a linear relation between the isoscaling parameter a, and the difference 
in the isospin asymmetry (Z/A) 2 of the fragments, with appreciably different slopes, 
was predicted for two different forms of the density dependence of the symmetry 
energy ; a " stiff " dependence (obtained from Gogny-AS interaction) and a " soft 
" dependence (obtained from Gogny interaction). 

Fig. 4 (left) shows a comparison between the experimentally observed a and 
those from the AMD model calculations plotted as a function of the difference in 
the fragment asymmetry for the beam energy of 35 MeV/nucleon. The solid and 
the dotted lines are the AMD predictions using the " soft " (Gogny) and the " 
stiff " (Gogny-AS) density dependence of the symmetry energy, respectively. The 
solid and the hollow symbols (squares, stars, triangles and circles) are the results 
of the present study for the two different values of the fragment asymmetry, as- 
suming Gogny and Gogny-AS interactions, respectively. Also shown in the figure 
are the scaling parameters (asterisks, crosses, diamond and inverted triangle) taken 
from various other works in the literature. It is observed that the experimentally 
determined a parameter increase linearly with increasing difference in the asym- 
metry of the two systems as predicted by the AMD calculation. Also, the data 
points are in closer agreement with those predicted by the Gogny-AS interaction 
(dotted line) than those from the usual Gogny force (solid line). The slightly lower 
values of the symbols from the present measurements with respect to the Gogny- 
AS values (dotted line) could be due to the small secondary de-excitation effect of 
the fragments not accounted for in this comparison. It has been shown [ 111] that 
the experimentally determined a values can be lower by about 10 - 15 % for the 
systems and the energy studied here. Accounting for this effect results in a slight 
increase in the a values bringing them even closer to the dotted line. The observed 
agreement of the experimental data with the Gogny-AS type of interaction therefore 
appears to suggest a stiffer density dependence of the symmetry energy. However, 
as mentioned in section 2, the effect of secondary de-excitation in the dynamical 
model calculations is currently under study [ I2()j and the predicted sensitivity may 
be significantly diminished by the secondary decay. 

Recently, Chen et al. [ I21| also showed, using the isospin dependent Boltzmann- 
Uehling-Uhlenbeck (IBUU04) transport model calculation, that a stiff density de- 
pendence of the symmetry energy parameterized as E sym w 31.6 {p/ p a ) lm explains 
well the isospin diffusion data [ 22 from NSCL-MSU (National Superconducting Cy- 
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Fig. 4. (left) Scaling parameter as a function of the difference in the fragment 
asymmetry for 35 MeV/nucleon. (Right) Parameterization of the various forms of 
the density dependence of the nuclear symmetry energy used in the analysis. 

clotron Laboratory at Michigan State University). Their calculation was also based 
on a momentum-dependent Gogny effective interaction. However, the present mea- 
surements on isoscaling gives a slightly softer density dependence of the symmetry 
energy at higher densities than those obtained by Chen et al. 

The difference in stifness is clear from figure 4 (right), which shows the pa- 
rameterization of various theoretical predictions of the density dependence of the 
nuclear symmetry energy in isospin asymmetric nuclear matter. The dot-dashed, 
dotted and the dashed curve corresponds to those from the momentum dependent 
Gogny interactions used by Chen et al. to explain the isospin diffusion data. These 
are given as, E ss/m 31.6 (p/po) 1 , where, 7 = 1.6, 1.05 and 0.69, respectively. The 
solid curves and the solid points correspond to those from the Gogny and Gogny-AS 
interactions used to compare with the present isoscaling data. As shown by Chen 
et al, the dependence parameterized by E sym w 31.6 (p/ Po) 1 ' 05 (dotted curve) ex- 
plains the NSCL-MSU data on isospin diffusion quite well. On the other hand, 
the isoscaling data from the present work can be explained well by the Gogny-AS 
interaction (solid points). 

Both measurements yield similar results at low densities with significant differ- 
ence at higher densities. It is interesting to note that by parameterizing the density 
dependence of the symmetry energy that explains the present isoscaling data, one 
gets, E sym « 31.6 (p/po) 1 , where 7 = 0.69. This form of the density dependence 
of the symmetry energy is consistent with the parameterization adopted by Heisel- 
berg and Hjorth- Jensen in their studies on neutron stars [EH- By fitting earlier 
predictions of the variational calculations by Akmal et al. [ I24| , where the many- 
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body and special relativistic corrections are progressively incorporated, Hcisclbcrg 
and Hjorth- Jensen obtained a value of E sym (p ) = 32 MeV and 7 = 0.6, similar 
to those obtained from the present measurements. The present form of the den- 
sity dependence is also consistent with the findings of Khoa et al. [ I25| , where a 
comparison of the experimental cross-sections in a charge- exchange reaction with 
the Hartree-Fock calculation using the CDM3Y6 interaction [EH], reproduces well 
the empirical half-density point of the symmetry energy obtained from the present 
work (see fig. 2 of Ref. [E>]). 

The observed difference in the form of the density dependence of the symmetry 
energy between the present measurement and those obtained by Chen et al. is not 
surprising. Both measurements probe the low density part of the symmetry energy 
and are thus less sensitive to the high density region. But the important point 
to be noted is that both measurements clearly favor a stiff density dependence 
of the symmetry energy at higher densities, ruling out the very " stiff " (dot- 
dashed curve) and very " soft " (solid curve) predictions. These results can thus 
be used to constrain the form of the density dependence of the symmetry energy at 
supranormal densities relevant for the neutron star studies. 

In view of the findings from the present measurements and those of Chen et 
al., we believe that the best estimate of the density dependence of the symmetry 
energy that can be presently extracted from heavy ion reaction studies is, E SJ(m 
31.6 (p/po) 7 , where 7 = 0.6 - 1.05. Measurements at higher densities should be 
able to constrain the density dependence of the symmetry energy further. 

5. Conclusions 

In conclusion, a number of studies have been carried out at TAMU to study the 
symmetry energy in the equation of state of isospin asymmetric nuclear matter. 
The results were analyzed within the framework of statistical and dynamical model 
calculations. It is observed that the properties of nuclei at excitation energy, isospin 
and density away from the normal ground state nuclei are significantly different and 
sensitive to the symmetry energy. The symmetry energy required to explain the 
isoscaling parameter of the fragments produced in multifragmentation reactions are 
significantly lower, and as small as 15 MeV. The dynamical model calculation of 
the isoscaling parameter shows that a stiffer form of the density dependence of the 
symmetry energy is preferred over a soft dependence. A dependence of the form 
E sym ~ 31.6 {p/po) ' 69 appears to agree better with the present data. Recently 
it has been shown that this form of the density dependence of symmetry energy 
provides an accurate description of several collective modes having different neutron- 
to-proton ratios. Among the predictions from this dependence are a symmetric 
nuclear-matter incompressibility of K = 230 MeV and a neutron skin thickness in 
208 Pb of 0.21 fm. Further, this dependence leads to a neutron star mass of M max 
— 1.72 Mq and a radius of R = 12.66 km for a "canonical" M = 1.4 M Q neutron 
star. These results have significant implications for nuclear astrophysics and future 
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experiments probing the properties of nuclei using beams of neutron-rich nuclei. 
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